Absence of surface mode in a visco-elastic material with surface tension 



00 



Hiizu Nakanishi and Satoshi Kubota* 
Department of Physics, Kyushu University 33, Fukuoka 812-8581, Japan 



The surface waves in the visco-elastic media with the surface tension are studied using the Voigt- 
Kelvin model of the visco-elasticity. It is shown that the surface mode of oscillation does not exist 
in the parameter region where the effect of surface tension is larger than that of the elastic stress at 
the surface unless the viscous stress masks the elastic stress in the bulk. In the region, the surface 
oscillation is suppressed and the oscillation beneath the surface diffuses after the pulse goes into the 
bulk. The experimental relevance of the present results is also discussed. 



3 

On 



o 



i 

C 

o 
o 



o 

00 



T3 

C 

O 
o 



Surface modes of oscillation in a material are quite 
different from bulk modes due to the boundary condi- 
tion at the surface. In an elastic material, the longitudi- 
nal and the transverse modes are mixed in the Rayleigh 
mode (surface mode), and its speed is slower than both 
of them. As for the fluid where only restoring force is 
the surface tension, the dispersion of the surface wave is 
not linear but given by u> oc fc 3 / 2 . The transition between 
the two modes in the visco-elastic material like polymer 
solutions or gels has been studied theoretically and 
experimentally |^-^]. 

The mode characters of these two modes, however, are 
quite different and they cannot mix easily in a visco- 
elastic material where both the elastic and the surface 
tension operate. It has been expected that the spectral 
peak of the surface tension wave in the thermal fluctua- 
tion does not transform into the Rayleigh peak smoothly 
but the two peaks tend to co-exist in the transition re- 
gion, and the spectral structure does not be simple 
The expected shoulder structure of the spectral peak has 
not been observed clearly but the substantial increase of 
the increase of the spectral width was demonstrated in 
the transition region for polymer solutions upon increas- 
ing density §,0 

Recently, a careful experiment was done to measure 
the speed and the attenuation of the externally excited 
surface wave in the sol and the gel phase of tungstic acid 
pof , which undergoes a sol-gel transition upon decreas- 
ing pH, and it was shown that there is an anomalous 
behavior around the sol-gel transition point; the wave 
speed rises sharply and then drops discontinuously, and 
the attenuation rate has a sharp peak. These behaviors 
have not been expected by the previous theories. 

Being motivated by the above experiment, we analyze 
the dispersion relation of surface modes of a visco-elastic 
material. We focus our attention to the experimental 
situation where the inverse of the surface wave frequency 
is much smaller than the structural relaxation time of 
the polymer network over the wave length. Then the 
Maxwell type of visco-elasticity can be ignored and the 
Voigt-Kelvin visco-elasticity, in which the stress is the 
sum of the elastic and the viscous stresses, is a good 
model. Note that the Voigt-Kelvin model does not be- 



come particularly bad around the sol-gel transition point 
because we are interested in the phenomenon with finite 
frequency and wave length. The sol phase behavior is 
also obtained within this model by setting the elasticity 
zero. 

With this model, we demonstrate that the surface 
mode does not exist around the region where the sur- 
face tension competes with the elastic stress. 

We assume the visco-elastic material is occupied in the 
z > space and the surface is located at the x — y plane. 
The displacement field u(x, y, z, t) can be decomposed 
into the transverse part u\ and the longitudinal part u*i 
as 



u = u t + ui; 



V • ut = 0, VxS, = 0, 



(1) 



then the equations of motion for u t and ui can be ex- 
pressed as 



/"'/ = [E + rj^] V-V, . 



pui 
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V 2 ui 
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where the dot denotes the time derivative and p is the 
mass density, E and K are the shear and bulk modulus, 
and r\ and £ are the shear and bulk viscosity, respectively. 
Using a potential function <pi, the longitudinal displace- 
ment ui can be expressed as 



then equation of motion (|^) is expressed as 
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The stress tensor cr^ is given by 
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Tr[e] <5, 
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where e is the strain tensor 



1 / dui duj 

2 \dxj dxi 



(7) 



For simplicity, we take the incompressible limit, K — » 
oo, which is good for the experimental situation where 
the surface wave speed under investigation is much slower 
than the longitudinal wave speed. Then we have 







and 



"n = -> | E + r i§i ) : " 



n Oij. 



(8) 



(9) 



Now we examine the solution of the form 

u t (x,z,t) = (-f(z,t)/ik, 0, /(M))e ifc *, (10) 

0,(1, z,t) =^)e 4fc ^ fc2 , (11) 

assuming the physical quantities does not depend on the 
y co-ordinate and the y component of the displacement 
is zero. In (|l0|), the prime denotes the derivative in z. 
The parameter k is the wave number in the x direction; 
k > in order that the solution should be the surface 
mode. Then the equation of motion (0) is given by 



pf(z,t) = [E + r, 



'dt 



The boundary conditions at z = are 

f) 2 n 

&xz = 0, a zz 



dx 2 



(13) 



where 7 is the surface tension. These can be represented 
as 



- (f"(0,t) + k 2 f(Q,t))+2k 3 <f>(t) =0, 

2 ! E + v m) ( /,(0, t] + + p ^ t] = 



(14) 



7 fc 2 (/(0,t)-fc^(t)). (15) 
If we assume the sinusoidal time dependence, 

f(z, t) = foe-™-™*, </>(€) = ^e-' l "\ (16) 

with 

k = \Jk 2 — pu> 2 /n(ui), m( w ) = E — iwrj, (17) 



then, from (|14|) and (|15[), we obtain the equation to de- 
termine the dispersion relation, 



2/; 2 - ^.V - 1 P 4^ - ^Jk 2 - J^l. (18) 



/i(w) 



Note that the branch with the positive real part should be 
taken for the square roots in order that the mode should 
be localized around the surface. 

It is easy to examine that this equation gives the 
Rayleigh wave dispersion ui — cny Ejp ■ k with cr = 
0.9553, for the case of the elastic material (E ^ 0, 7 = 
r] = 0) and the capillary wave dispersion u> = y/^f/p • k 3 ^ 2 
in the case of the fluid (7 ^ 0, E = 77 = 0). 

In order to see how the Rayleigh mode and capillary 
mode exclude each other, we consider first the case with- 
out the dissipation (77 = 0, E ^ and 7 ^ 0), then 
eq.(jl^) becomes 



(2 



4VT 



(19) 



The dispersion relation is obtained by determining c for a 
given value of k. Again, note that the square root should 
have a positive real part. 

This equation gives the Rayleigh wave like dispersion 
for k < E/j, In the case of k > E/j, however, it does 
not have a solution in the physical branch of the square 
root. This means that there is no surface mode in the 
region of k > E/j, where the capillary stress dominates. 

In order to examine what is happening in the region 
without the surface modes, we solve the equation of mo- 
tion ([12]) numerically with 7/ = under the boundary 
conditions ( |l4| ) and ( [l5|) from the initial state given by 
f(z,t = 0) localized around the surface within a length 
A as, 



/(z,t = 0) = / o e- z/A . 



(20) 



The time developments are shown in Fig.l, where energy 
density, 



1 / duj du q 
-E 1 3 



Z,J=X,Z 



4 V dxj dxi 



(21) 



is plotted as a function of z for k = 0.5, 1, and 5 (E/ r y) 
with the initial configuration with A = 1, 0.5, and 0.5 
(7/-E), respectively. Note that the initial shape of u e 
is not proportional to (e~ z ^) 2 because the longitudinal 
component given by (H) and ( O ) decays as e~ kz in z. 

It can be seen that the surface oscillation remains af- 
ter the pulse is gone into the bulk for the case of k = 
0.5 (E/j), where the surface mode is allowed. On the 
other hand, the oscillation at the surface is strongly sup- 
pressed for k = 5 (.E/7), where the surface mode is not 
allowed. It is interesting to see, in this case, that there is 
a peak of oscillation just beneath the surface in addition 
to the pulse that travels into the bulk, and the peak of 
the oscillation diffuses away as the pulse goes in the bulk. 
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FIG. 1. Time developments of the oscillation for the rj — case. The energy density is plotted against z in the unit of 
(■y/E) for the wave numbers in the x direction k =0.5, 1, and 5 (E/y) with the initial state given by (GOJ) with A = 1, 0.5, and 
0.5 (y/E), respectively. The time is shown in unit of (p/E) 1,/2 y/E. 



If the viscosity is introduced into the system, the sit- 
uation becomes slightly more complicated. There takes 
place the competition among three stresses; the elastic 
stress, the surface tension, and the viscous stress. De- 
pending upon which stress dominates, the mode charac- 
ter is expected to be one of the followings, the Rayleigh 
mode, the capillary mode, or the overdamped liquid 
mode. In addition to these, however, there is the pa- 
rameter region where no surface mode is allowed, as we 
have seen. 
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FIG. 2. Dispersion relations for a = E(r\/y) 2 / p = 0.5, 
0.3, and 0. The real part of Q. = (rf / py 2 )u> is plotted against 

Q = (v 2 / P /y)k. 

In Fig. 2, the result of dispersion relation for E = 
0.5, 0.3, and pid/f]) 1 is given; the real part of £1 = 
(if / frf 2 )u is plotted as a function of Q = (rj 2 / fq)k. In 
the case of E = 0.3 pij/il) 2 , for small Q, the dispersion 
is linear and the mode is the Rayleigh wave like, but for 
an intermediate value of Q, the mode becomes capillary 
wave like, and for large Q, £1 becomes pure imaginary 
and the mode is the overdamped liquid mode. Between 
the elastic mode and the capillary mode, there exists the 



region where there exists no surface mode, namely, there 
is no solution for Q on the physical branch of the square 
root in eq. (jl^) ; we call it the "gap region" . 

The phase diagram in the E — k plane has been ob- 
tained numerically by calculating the dispersion from 
( |l8| ) [ fPH , but the result is illustrated in Fig. 3 schemat- 
ically. It can be understood as in the following way. 
The boundary between the Rayleigh mode and the over- 
damped liquid mode is determined by the competition 
between the elastic stress and the viscous stress, namely, 
Ek\u\ and r]Luk\u\, therefore the boundary is given by 

k ~ -JpE, (22) 
V 

because u> ~ y/ E / p ■ k in the Rayleigh mode region. Sim- 
ilarly, the boundary between the capillary wave and the 
overdamped liquid mode is given by 

k~pj/ V 2 . (23) 
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FIG. 3. Schematic phase diagram in the E — k plane in 
the logarithmic scale. The shaded region is the gap region, 
where the surface wave does not exist. 
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It is instructive to examine how the gap region arises 
between the Rayleigh and the capillary region. For small 
enough k, the elastic stress dominates and the surface 
wave is Rayleigh like. For E < ph/i]) 2 , upon increasing 
k, the surface tension becomes dominant at 



k ~ Eh, 



(24) 



where the surface mode disappears as we have seen for 
the dissipationless case. As k increases further, the elas- 
tic stress is taken over by the viscous stress at 



p E 
7 V 



2/3 



(25) 



then the capillary wave emerges 12 . 

It is very suggestive that the existence of the capillary 
wave depends upon the competition between the elas- 
tic stress and the viscous stress, not the surface tension. 
No matter how large the surface tension is, the capil- 
lary wave is not allowed under the existence of elastic 
stress unless there exists the viscous stress. This can be 
explained as in the following: The mode characters of 
the two are too different to form a mixed mode, there- 
fore, the Rayleigh wave dominates rather than mixes with 
the capillary mode when the elastic stress is larger. On 
the other hand, the capillary wave cannot overtake the 
Rayleigh wave even when the surface tension at surface is 
larger than the elastic stress because the surface tension 
operates only at the surface and the displacement field 
inside is governed by the elastic stress in the bulk. The 
viscosity, however, gives the capillary wave a chance to 
win by masking the effects of the elastic stress. 

The overdamped liquid mode under the existence of 
the surface tension also deserves some comments: In the 
bulk, it has the dispersion 



^-k\ 
P 



(26) 



As for the surface mode, there is a similar overdamped 
mode with the dispersion 



to 



-% - c R k 
P 



(27) 



with the same constant cr with the Rayleigh wave, and 
this is the only one in the case without the surface ten- 
sion. With the surface tension, however, there appears 
another mode, whose damping rate is smaller than the 
other's; the dispersion is given by 



CO 



-i — k 
277 



(28) 



for k 3> ip/rj 2 . 

Let us examine the experimental relevance of this gap 
region of the surface wave. In the recent experiment |L0| ], 
it was reported that the surface wave shows anomalous 
behavior in the tungstic acid; sharp increase and discon- 
tinuous drop of wave speed and attenuation peak of the 



surface wave in the frequency range ~ 100Hz have been 
observed around the sol-gel transition point, where the 
elasticity is very small and the Rayleigh wave speed is 
of order of 0.1 m/s. Within the present framework, the 
sharp increase of the wave speed can be interpreted as 
the result of the group speed increase at the upper edge 
of the gap in the dispersion curve. The discontinuous 
drop of the wave speed and the large attenuation of the 
surface wave come from the suppression of the surface 
oscillation in the gap region. 

The typical values of the parameters for the material 
are p ~ 10 3 kg/m 3 , r) ~ l(T 2 Pa • s, 7 ~ lO^N/m. The 
elastic modulus changes rather rapidly around the sol-gel 
transition point, but E ~ lOPa in the region where the 
sound speed is 0.1m/ s, then the expected gap region is 



or lu 



10 1 



This is 



consistent with the region where the anomaly observed. 

Before concluding, let us discuss possible effects of 
structural relaxation, which we ignore in the present 
work. The structural relaxation will relax part of the 
elastic stress for the longer time scale than its relaxation 
time r s t, and can be expressed by the Maxwell model of 
visco-elasticity in a phenomenological level. If we include 
the effect, the capillary mode region would extend into 
the gap region to some degree in the small k side because 
the structural relaxation reduces the bulk elastic stress in 
longer time scale than the structural relaxation time r s t • 
The way how the gap region shrinks from the present re- 
sult (Fig. 3) depends upon detailed structure of the model 
and the relevant time scale; for to 3> 1/T st (lj <C l/r s t), 
the part of the stress can be treated as elastic (viscous). 
The basic feature of the phase diagram, however, should 
not change as long as there exists a pure elastic compo- 
nent of the stress in the system as we assumed in the 
present study. 

In summary, we found that, in the visco-elastic mate- 
rial, there is a parameter region where the surface mode 
is absent and it corresponds to the region where the sur- 
face wave anomaly is found in the recent experiment on 
tungstic acid |l(J. 

The authors would like to thank H.Okabe, 
K.Kuboyama, K.Hara, and S.Kai for showing their ex- 
perimental data prior to publication. This work is par- 
tially supported by Grant-in- Aid for Scientific Research 
(C) (# 0940468) provided by the Ministry of Education, 
Science, Sports, and Culture, Japan. 



* Present address: NTT Network Service Systems Labora- 
tories, 3-9-11, Midori-Cho, Musashino, Tokyo 180, Japan 
[1] H. Pleiner, J.L. Harden, and P. Pincus, Europhys. Lett. 
7, 383 (1988). 

[2] J.L. Harden, H. Pleiner, and P.A. Pincus, J. Chem. Phys. 
94, 5208 (1991). 



4 



[3] J. Jackie and K. Kawasaki, J. Phys.iCondens. Matter 7, 
4351 (1995). 

[4] H. Kikuchi, K. Sakai, and K. Takagi, Jpn. J. Appl. Phys. 

30, L1668 (1991). 
[5] H. Kikuchi, K. Sakai, and K. Takagi, Phys. Rev. B 49, 

3061 (1994). 

[6] R.B. Dorshow and L.A. Turkevich, Phys. Rev. Lett. 70, 
2439 (1993). 

[7] B.H. Cao, M.W. Kim, and H.Z. Cummins, J. Chem. 

Phys. 102, 9375 (1995). 
[8] H. Takahashi and P-K. Choi, Jpn. J. Appl. Phys. 35, 

2939 (1996). 



[9] K. Motonaga, H. Okabe, K. Hara, and K. Matsushige, 
Jpn. J. Appl. Phys. 33, 2905 (1994); ibid. 33, 3514 
(1994). 

[10] H. Okabe, K. Kuboyama, K. Hara, and S. Kai, Rept. 
Progr. Polym. Phys. Jpn. 40, 219 (1997); H. Okabe, K. 
Kuboyama, K. Hara, and S. Kai, Jpn. J. Appl. Phys. 37, 
2815 (1998). 

[11] S. Kubota and H. Nakanishi, Prog. Theor. Phys. (Kyoto) 

Suppl. 126, 359 (1997). 
[12] The gap region is missing in the phase diagram Fig.l of 

Q| , where the boundary between the capillary mode and 

the elastic mode is drawn. 



5 



